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ABSTRACT 

Over cosmic time, galaxies grow through the hierarchical merging of smaller galaxies. However, 
the bright region of the galaxy luminosity function is incompatible with the simplest version of hi- 
erarchical merging, and it is believed that feedback from the central black hole in the host galaxies 
reduces the number of bright galaxies and regulates the co-evolution of black hole and host galaxy. 
Numerous simulations of galaxy evolution have attempted to include the physical effects of such feed- 
back with a resolution usually exceeding a kiloparsec. However, interactions between jets and the 
interstellar medium involve processes occurring on less than kiloparsec scales. In order to further the 
understanding of processes occurring on such scales, we present a suite of simulations of relativistic 
jets interacting with a fractal two-phase interstellar medium with a resolution of two parsecs and a 
largest scale of one kiloparsec. The transfer of energy and momentum to the interstellar medium 
is considerable, and we find that jets with powers in the range of 10 43 -10 46 erg s _1 can inhibit star 
formation through the dispersal of dense gas in the galaxy core. We determine the effectiveness of 
this process as a function of the ratio of the jet power to the Eddington luminosity of the black hole, 
the pressure of the interstellar medium and the porosity of the dense gas. 

Subject headings: galaxies: evolution - galaxies: formation - galaxies: jets - hydrodynamics - ISM: 
jets and outflows - methods: numerical 



1. INTRODUCTION 

It is widely believed that feedback from active galactic 
nuclei (AGN) during the epoch of galaxy formation is 
required to explain the relation between black hole and 
bulge mass/velocity dispersion ([M agorria n^ et al.| |1998| 
IGebhardt et al.||2QQQ| |Tremaine et al.||2002| ), the d ericit 
ol bri ght gala xiesin tne galaxy luminosity function ( |Cole| 



et al.||2001| |Norberg et aT|2002| |Huang et alT]|2QQ3[ )rand 

the completion ol star formation m massive galaxies at 
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the infall of star-forming gas once t he central black hol e 
grows to a cri t ical size. Ac cordingly, Silk & Rees (1998) 



Fabian (1999), King (2005), and others, have appealed to 
the physics of either energy-driven or momentum-driven 
bubbles in order to explain the relationships between the 
mass of the black hole and the parameters of the host 
galaxy. 

In order to m odel the galaxy luminosity function, Cro- 
ton et al.| ( |2006[ ) have utilized semi-analytic models based 



on the output of the Millennium Simulation, incorporat 
ing "radio-mode" feedback, coupled with a prescription 
for the accretion rate into the center of each evolving 
galaxy. Their feedback prescription is motivated by the 
well-documented evidence for the effect of radio galaxies 
on "coolin g flow" galaxies (Fabian et al. 2003 McNamara 
et al.||2005| ). 



There is also a growing literature on cosmological sim- 
ulations in a ACDM cosmogony involving both dark 
matter and gas dynamics, which incorporate feedback 
from both supernovae (SN) and black holes, a nd which 
test scenarios of galaxy merging and growth ( Springel 
fc Hernquistl|2003b|a| |Booth fc Schaye||2009b|a| |Schaye 



et ah1|20 10|. In simulations using the Gadget-2 code 



(Springel 2005) the total number of SPH particles repre- 
senting the baryonic component exceeds 250 million, and 
the effective spatial dynamic range is an equally impres- 
sive 10 5 per dimension. Nevertheless, the best spatial 
resolution is about 2 kpc and does not resolve the spa- 
tial scales where important dynamical processes occur. 
This is highlighted by the prescriptions for black hole 
gr owth and feedback describ ed, for example in the work 
by Booth fc Schaye| (2009a). Accretion is described in 
terms of the Bondi-Hoyle accretion rate multiplied by a 
factor, which can be as large as 10 2 . The rationale for 
this approach is that at sub-grid scales the density would 
be larger, and the real accretion rate would be apprecia- 
bly higher. However, the higher densities and the conse- 
quent cooling and fragmentation on sub-grid kpc scales 
would create a multi-phase interstellar medium (ISM), 
the physics of which is not satisfactorily captured by 
the simulations. In particular, such a medium is porous, 
and the interaction between jets and the dense, poten- 
tially star-forming clouds of gas is complex, with radio- 
emitting plasma being able to channel through holes in 
the density distribution, rather than isotropically im- 
pacting a smo oth distribution of dense gas as shown in 
previous work ( [Sutherland fc Bicknell|[2007l ). 

Given that black- hole driven feedback occurs in bright 
galaxies, there does not appear to be a consensus on the 
type of feedback: How much AGN power is involved, 
and does the feedback involve radiative or me chanical 
proce s ses or both? In their radio-mode model, |Croton 



et al. ( |2006[ ) attribute feedback to radio galaxies accret- 
ing at rates well below Eddington, and for typical el- 
lipticals, this means low-powered Fanaroff- Riley Class I 
radio galaxie s are the primary drivers of feedback. In 
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Fabian (1999]) the momentum for disper sing the cir cum- 
nuclear gas comes from a quasar wind; in King (2005 ) the 



2 



Wagner & Bicknell 



momentum of the outflow is provided by an Eddington- 
limited radiation-driven wind; in SPH simulations the 
accretion rate, which is the ultimate power s ource for 



scalars in the relativi stic hydro dynamic solver. Details 



an outflow, c an approach Eddington values (Booth & 
|Schaye||2009b| ). 

There is also an issue of what class of AGN actually 
drives black hole feedback. The models by [Croton et al.| 
(2006) invoke low-powered radio galaxies! Observers, 



however, have focused on powerful radio galaxies, mainly 
at z > 2; in these galaxies there is evidence for substan- 
tial outflows of li ne-emitting gas and neutral gas driven 
by th e radio jets (Nesvadba et al. 2009; M organti et al.| 
2010). It is possible that there is a role lor radio galax- 



les with a range of powers: Powerful radio galaxies may 
be responsible for the establishment of the stellar mass 
at z > 2, and less powerful sources may be responsible 
for the maintenance of the stellar content th rough the 
inhibition of cooling flows ( jNulsen et al.|[2009[ ). We also 
note that radio galaxies are mainly relevant to the ellip- 
tical galaxy population, and that the separ ate luminosity 



functions for early and late- type galaxies (Huang et al. 
2003 ) indicate the requirement for feedback in both pop- 
ulations. 

In this paper we consider the potential role of radio 
galaxies in AGN feedback and address the following ques- 
tions: (1) What jet power is required for the radio galaxy 
phase to have an important effect on inhibiting star for- 
mation in a given host, and (2) Is the range of radio 
powers broad enough that radio galaxies could affect the 
entire distribution of bright ellipticals? In this paper we 
present progress in answering these questions through 
simulations with a resolution of 2pc per computational 
cell. These simulations confront the sub-grid physics that 
current large scale SPH simulations do not address. In 
particular, we consider the effect of powerful relativistic 
jets on a two-phase ISM consisting of hot gas, in which 
is embedded a dense porous phase of warm gas. 

These simulations extend the simulations described in 



(Sutherland & Bicknell 2007), in which a jet with a ki- 
netic power of 3 x 10 43 ergs _1 propagates through an 
inhomogeneous medium in the form of an almost Kep- 
lerian fractal disk. It is evident from that simulation 
that in the geometry considered, jets of that power could 
not exert enough impact on the clouds to disperse them, 
and that the jets would not have a important effect on 
star formation in the core of the host galaxy. In the 
present simulations we consider jets with powers rang- 
ing from 10 43 to 10 46 ergs -1 propagating through a two- 
phase medium, in which the dense clouds are spherically 
distributed throughout a region 1 kpc in diameter. These 
initial data are meant to describe a typical protogalaxy, 
in which dense gas has accumulated in the core. The 
fractal distribution of the dense gas enables us to di- 
rectly examine the effect of porosity on the evolution of 
potentially star-forming clouds. 

In the following sections we describe the parameters 
of the simulations in more detail and then discuss our 
results. 

2. MODEL PARAMETERS AND INITIAL CONDITIONS 

In our simulations, we use the publicly available, open- 
source code FLASH (iFryxell et al.||2000[ ) version 3.2, 
to which we have added code to incorporate radiative 
cooling of thermal gas and code to advance advected 



% 

ties ot FLASH, utilizing up to seven levels of refinement; 
this corresponds to a nominal cubical simulation grid of 
1 kpc 3 in physical dimensions, consisting of 512 3 cells at 
maximum resolution. However, note that a restricted 
one parameter scaling of phy sical dimensions is possible 
( [Sutherland fc Bicknell][2007| . 

The relevant jet parameters, which are initiated and 
maintained constant at the boundary that is the jet 
inlet, are: the Lorentz factor Fj et = (1 — /3? et ) -1 / 2 , 
where the /3j et = v- ]et /c is the jet velocity in units of 
the speed of light; the proper density parameter, x — 
(7 - 1) PjetC 2 /7Pjet, where 7 is the polytropic index, p jet 
is the rest mass density, and p- ]et is the pressure; and the 
jet power. Let Aj e t be the jet cross-sectional area. The 
jet power is 



of the solve r are in (Mignone & Bodo 2005 Mignone 

"11 

of 



et al.l 20051) . We exploit the adaptive mesh capabili- 
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jet = ^— yCp jet r jet /3 je tA jet 
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(i) 



We adopt 7 = 5/3 for both jet and ambient gas. In 
all simulations presented here, Fj et = 10, and x — 1-6- 
The jet inlet is a circular region of area Aj et centered 
at (x,y,z) = (0,0,0) with normal (1,0,0). The ini- 
tial jet velocity is parallel to the x-axis. The boundary 
x = is reflective, apart from the jet inlet. All other 
boundaries of the simulation domain are designated as 
inflow / outflow boundaries . 

The two parameters describing the hot phase of the 
ISM are the temperature, X^, which is fixed at 10 7 K, 
and the total number density n^. The pressure, p/k, 
where k is Boltzmann's constant, for the ISM in each of 
the simulations is either 1 6 or 10 7 , typical of the ISM 



in giant elliptical galaxies ( Mathews fc Brighenti| |2003 ) . 
The distribution of den se, warm clouds is prescribed m 



similar fashion to that in Sutherland fc Bicknell| ([2 007), 
follow ing work on terrestrial clouds by |Lewis fc A ustin 
(2002). The main difference is that, here, the average 



mass dens ity is uniform, whereas in [Sutherland fc Bick- 
|nell| ( |2007| ) the average density is that of a near- Kepler ian 
disk. For completeness, we describe the other details of 
the cloud distribution used in this series of simulations: 

(i) We begin by constructing a cube of density fluc- 
tuations in which the single point statistics of the 
density distribution are described by a log-normal 
distribution. The mean, /i, of the parent distribu- 
tion is 1.0, and the variance, a 2 = 5.0, these values 
being consistent with starburst reddening and ex- 



tincti on models ( Fischera et al.|2003| |Dopita et al. 
2004). This valued" a 1 is the same as that used 



in Sutherland & Bicknell 



( |2QQ7] ). 

:bulenc( 



In numerical sim- 



ulation s ol supersonic turbulence, |Federrath et al. 



(2010) find a 2 ~ 3.6 and 35 for solenoidal (diver- 



gence free) and compressive (curl free) forcing re- 
spectively, so that our adopted value of 5 is closer 
to their solenoidal result. 

(ii) The power spectrum of the density distribution 
D(k) oc k~ b ^ for the range of waves numbers 



< k < k n 



and zero outside of this range. 



The parameter k m i n sets the maximum cloud size. 



Relativistic Jet Feedback 



3 



For these simulations k m [ n = 20 (in cell units in 
Fourier space) and this value limits the maximum 
size of an individual cloud to approximately 25 pc 
for 512 cells along each axis of the lkpc compu- 
tational cube. Limiting the cloud size to 25 pc al- 
lows for appropriate coverage and variation of hot 
phase ISM and warm phase ISM along any path 
through the region of our simulation grid filled with 
clouds. This ensures that the spatial distribution of 
clouds is approximately isotropic, and that the jet 
plasma encounters a statistically significant num- 
bers of clouds along its main axis of propagation. 
As a result, the flow of the jet plasma through the 
hot phase is approximately isotropic. The parame- 
ter /c max is set by the resolution of the simulation. 
Here & max = 255 and is equivalent to two compu- 
tational cells. 

(iii) The parent, unit mean distribution is scaled by the 
mean particle density of warm clouds, (n w ); this 
parameter is determined by the ratio of (n w ) to 
n h - 

(iv) The temperature in each cell of the warm gas distri- 
bution is determined by pressure equilibrium with 
the hot gas. When the warm gas temperature ex- 
ceeds T cr it = 3 x 10 4 K (in the lowest density parts 
of the distribution), it is deemed to be thermally 
unstable and is replaced by hot gas. This makes the 
gas porous. We can then define a volume filling fac- 
tor, /y, of the warm phase by integrating over the 
log-normal density distribution above the critical 
density p crit = p upi$ M / T CT - lt k , where ft w 0.6156 
is the mean molecular weight and u is an atomic 
mass unit (se e also Appendix B of Sutherland & 
BicknelT|[2007| ). 



For a given value of T crit , the volume filling factor 
is determined by the ratio of the mean warm phase 
density to hot phase density, (n w ) juh\ the larger 
this parameter the larger the filling factor. This 
log-normal, fractal, porous distribution of clouds is 
an important feature of the simulations presented 
here. 



(2009) identi- 



Silk|fll997}|200lD and |Silk fc Norman 
tied porosity as a key parameter in their models of 
feedback-regulated star-formation, since the poros- 
ity determines the extent to which energy-driven 
bubbles generated by SN or AGN activity are con- 
fined. The porosity in the two-phase ISM employed 
in our simulations plays a similar role, in that it 
constrains the progress of the jet along its principal 
axis and ensures the confinement of the pressurized 
bubble for times much larger than the dynamical 
time of an unimpeded jet. 

The porosity of dense clouds is the primary rea- 
son for the differences in the evolution of the radio 
source and the differences in energy and momen- 
tum imparted by the jet to the ISM between our 
simulations and those of a jet propagating into a 
uniform medium. 



(v) The radial extent of the gas is truncated by a sphere 
of radius 0.5 kpc. Hence, we are simulating the 



interaction of jets with warm gas within 1-2 core 
radii of the parent galaxy. This justifies the neglect 
of gravity. 

(vi) We include non-equili brium, optically thin atomic 
cooling for T > 10 4 K fl Sutherland fc Dopita| l993" 



Sutherland et al. 2003); for T < 10 4 K the cooling 



is set to zero. We u se updated solar abundances 
( |Asplund et"aL|2005| >. 



The non-equilibrium, cooling function is pre- 
calculated for a high temperature (10 10 K) shock. 
The main reason for using such a function was out- 
lined in § 2.1 of ( [Sutherland et al.||2003[ ) and, for 
completeness, we summarize it here: Following a 
shock there is a short-lived cooling spike as low- 
ionization material is suddenly shocked to a high 
temperature. In a numerically well-resolved shock 
this phase is short-lived and provides about 1% of 
the total cooling. However, in a simulation in which 
each shock is not well-resolved (as is the case here) 
this spike can dominate the cooling and the situa- 
tion is exacerbated by the interpolation of interme- 
diate temperatures at each cell. Apart from the ini- 
tial (and unimportant) cooling spike, our adopted 
non-equilibrium cooling function is a good approx- 
imation to the cooling function of shocks for veloc- 
ities above 150kms _1 . Since the cooling function 
is pre-calculated and interpolated within the code 
there are no additional computational costs associ- 
ated with a non-equilibrium calculation. 

The utilization of an optically thin cooling function also 
requires some comment, given the dense, reasonably ex- 
tended regions that evolve within the simulations. Con- 
sider a region with electron density n e = 10 3 n Gj s cm -3 
and thickness I. The electron scattering optical depth 
r e ~ 2.0 x 10 -3 n e ,3(//pc). In the simulations with the 
densest gas the initial average electron density n e ^ ~ 1 so 
that the largest clouds with I ~ 50 pc would be optically 
thin but verging on optically thick where the density is 
a factor of 10 above average. As a result of the radia- 
tive shocks, the electron density increases by a factor of 
10 — 100. Hence, even some of the less dense regions could 
become optically thick to scattering. Scattering alone is 
not enough to invalidate the assumption of optically thin 
cooling, and an admixture of dust is required to cause 
absorption. Moreover, the cooling in such regions occurs 
before maximum density has been obtained, and photons 
emitted by cooling plasma following a radiative shock 
have an escape route through the lower density regions 
of the shock. In addition, the temperature of the gas 
in the cooling, optically thin region of a radiative shock 
would be so high that dust would be destroyed. Hence, 
optically thin cooling is a reasonable first approximation, 
which may be limited in the largest clouds. More refined 
radiation hydrodynamic simulations, possibly involving 
larger clouds, may need to take optical depth effects into 
account. 

In summary, the key model parameters used in the 
different simulations are the jet power and the densities 
of the hot and warm ISM phases, which determine the 
cloud filling factor. Table [l] summarizes the assigned and 
derived parameters that we used for our simulations. 
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Table 1 

Simulation parameters 



Simulation 


log Pjet (a) 

(erg) 


(cm" 3 ) 


PismA (c) 
(cm- 3 K) 


(cm" 3 ) 




M™,tot (f) 
(1O 9 M ) 


A 


45 


0.1 


10 6 








B 


46 


1.0 


10 7 


1000 


0.42 


16 


C 


46 


0.1 


10 6 


100 


0.42 


1.6 


C 


46 


0.1 


10 6 


30 


0.13 


0.32 


D 


45 


1.0 


10 7 


1000 


0.42 


16 


D' 


45 


1.0 


10 7 


300 


0.13 


3.2 


E 


45 


0.1 


10 6 


100 


0.42 


1.6 


E' 


45 


0.1 


10 6 


30 


0.13 


0.32 


F 


44 


0.1 


10 6 


100 


0.42 


1.6 


F' 


44 


0.1 


10 6 


30 


0.13 


0.32 


G 


44 


1.0 


10 7 


1000 


0.42 


16 


G' 


44 


1.0 


10 7 


300 


0.13 


3.2 


H 


43 


0.1 


10 6 


100 


0.42 


1.6 



Note. — Runs labeled with primed ( u/ ") letters denote lower filling factor 
counterparts to runs with the same letter. 

( a ) Jet power 

( b ) Density of hot phase 

( c ) p/k of both hot and warm phases 

( d ) Average density of warm phase 

( e ) Volume filling factor of warm phase 

( f ) Total mass in warm phase 

3. RESULTS 

3.1. Morphology of radio source and interstellar medium 

The panels in Figs, [I] and [2] show the evolution of the 
density in the two simulations A and E, respectively. The 
figures present slices through the mid-plane z = of the 
grid. A is a reference simulation of a jet and homo- 
geneous ISM. The progress of such a jet has been well 
established over two decades of research, although we do 
note the unstable jittering of the t hree-dimensional jet 
once it has tra versed about 0.5 kpc (fMizuno et al.| |2007, 
also evident in [Sutherland fc Bicknell|2UU7| ), which leads 
to multiple hot spots. 

The progress of the jet in Fig. [2] (simulation E) is com- 
pletely different: Initially the jet is deflected in various 
directions as it floods through the porous screen of dense 
clouds, finding channels of least resistance and gradually 
dispersing the clouds through the effect of the ram pres- 
sure of the non-thermal plasma. In addition, the slower 
rate of progress of the jet traps the high pressure co- 
coon material; a pseudo-spherical bubble is driven into 
the ISM, and the dense clouds are driven outwards and 
dispersed as this bubble expands. Some of the dispersed 
material is accelerated to speeds ~ 1000 km s -1 while the 
speed of the densest parts of the clouds reaches several 
hundred kms -1 . 

The lateral extent of the cocoon of the A and E jets 
is very different. When the cloud- free jet has progressed 
about 0.4 kpc its radial extent is about 0.2 kpc, whereas 
the radial extent of the cocoon from the cloud-impeded 
jet is about 0.4 kpc. For the same length the latter jet 
processed a factor of 8 larger volume of the ISM, re- 
moving one of the reservations about the relevance of 
powerful radio galaxies to AGN feedback. Note also 
the streams of low density, pressurized cocoon material 
shooting out in random directions. This is a direct result 
of the porous ISM, into which the main jet is propagat- 
ing. Figures 2.1-2.5 and corresponding mpeg animations 



depicting the density evolution of runs E, C', D', G, and 
H, respectively, are available in the electronic edition of 
the Astrophysical Journal. 

We have generated images of synthetic radio surface 
brightness, based on the density of the non-thermal 
plasma and also assuming that the magnetic field is 
proportional to the pressure and a spectral index of 
0.6. A set of snapshots from simulation E are shown in 
Fig. [3j The interaction between jet and clouds produce 
a markedly different morphology from that of a classi- 
cal radio galaxy. Images such as these can be informa- 
tively compared with images of young GPS a nd CSS ra- 
dio galaxies (cf. Sutherland & Bicknell 2007). 



3.2. Velocity of dispersed clouds 

A convenient approach for assessing the effects of AGN 
feedback on galaxy formation is to examine the velocity 
imparted to gas, which could form new stars. The usual 
criterion for inhibition of further galaxy formation is that 
the velocities imparted to the clouds is gre ater than the 



veloc i ty dispersio n, a, of the host galaxy (|Silk fc Rees 



1998 King 2005). This does not necessarily mean that 



the clouds would be ejected to large radii but it does 
mean that the clouds would be highly dispersed within 
the potential well of the host galaxy. 

In order to track the different gas components we use 
tracers, which are the mass concentration of that com- 
ponent in each cell. In particular we use a warm gas 
tracer, <j) Wl which is initialized to unity in each cell of 
warm thermal gas. Figure [4] shows the velocity of the 
dispersed warm interstellar material, for <p w > 0.9, in 
simulation E at various epochs. At t ~ 98kyr the max- 
imum mean cloud velocity is reached. There is a wide 
range in velocity with the remaining dense cloud cores 
showing the lowest velocities, material resulting from 
the dispersal of the clouds showing intermediate veloc- 
ities ~ 300 kms -1 , and lighter material showing even 
higher velocities ~ 1000 kms -1 . Figures 4.1-4.5 and 
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Figure 1. Evolution of density in simulations A. The width and height of each panel are 1 kpc. This figure is also available as an mpeg 
animation in the electronic edition of the Astrophysical Journal. 




Figure 2. Evolution of density in simulation E. The width and height of each panel are 1 kpc. This figure is also available as an mpeg 
animation in the electronic edition of the Astrophysical Journal. Figures 2.1-2.5 and corresponding mpeg animations depicting the evolution 
of the density in runs E, C ; , D ; , G, and H, respectively, are also available in the electronic edition of the journal. 
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Figure 4. Speed of the warm gas for the value of the warm gas tracer <f> w > 0.9 in simulation E. The width and height of each panel 
are 1 kpc. This figure is also available as an mpeg animation in the electronic edition of the Astrophysical Journal. Figures 4.1-4.5 and 
corresponding mpeg animations depicting the evolution of the speed of the warm gas in runs E, C ; , D x , G, and H, respectively, are also 
available in the electronic edition of the journal. 



Relativistic Jet Feedback 



7 



corresponding mpeg animations that depict the veloc- 
ity evolution of the warm material, for which <p w > 0.9, 
in runs E, C, D', G, and H, respectively, are available in 
the electronic edition of the Astrophysical Journal. 

Taking p, the density in a cell, and v r , the radial veloc- 
ity in a cell, we compute for each simulation the mass- 
weighted mean radial velocity over all N cells defined 

by 

N 

Ys^w PVr 



(^V,u>) — 



1=1 



N 

E< 

1=1 



(2) 



We adopt (v r ^ w ) > a as the criterion for inhibition of fur- 
ther galaxy formation. This raises the question: What 
value of o~ is appropriate? If we use a value relevant to a 
giant elliptical a ~ 300 km s -1 , this discriminates against 
lower mass galaxies (a ~ 200kms _1 say), in which jet 
inhibition of star formation may be important in ear- 
lier epochs of galaxy formation. In order to establish the 
relevance of these simulations for all phases of galaxy for- 
mation we adopt the following procedure, in which the jet 
power is parametrized by its ratio 77 relative to the Ed- 
dington luminosity, that is, Pj e t = ^47tGMbh^pC<j^ 1 , 
where m p is the proton mass and cft is the Thomson 
scatt ering cross section. Adopting the Magorrian rela- 
tion ( |Tremaine et al.||2002| ) between black hole mass and 
velocity dispersion, Mbh ~ 8.1 x 10 6 af 00 M , where aioo 
is the velocity dispersion in units of 100 km s -1 , gives the 
jet power Pj et = 1.02 x 10 45 77 cr 4 00 ergs -1 , and the veloc- 
ity dispersion expressed in terms of the jet power in units 
of 10 45 erg s -1 is 



jet,45 1 



(3) 



We plot the results from each simulation on a dia- 
gram, Fig. [5j of the maximum value of (v r , w ) vs Pj et . 
Also plottea are the loci of a vs Pj et for incremental 
values of log 10 77, using equation ([3]). For a given value 
of 77 a jet will disperse the warm cloudy material if the 
corresponding ((v r ^ w ) ,Pj e t) point lies above the specific 
77 = constant locus. For example, consider the point for 
simulation F represented by the point Pj et = 10 44 erg s -1 , 
( v r,w) = 500kms _1 . If 77 = 10 -2 , the cloud speed 
comfortably exceeds the corresponding velocity disper- 
sion of approximately 180kms _1 . On the other hand, if 
rj = 10 -4 then the inferred velocity dispersion is approx- 
imately 560 km s -1 , and it is marginal whether feedback 
from such a jet would influence the evolution of the host 
galaxy. For 77 = 10 -5 , we expect no substantial feedback. 

A number of patterns are immediately obvious from 
this diagram. Consider first the sequence of simulations 
CEFH. These have the same values of Pism A ano ^ warm 
cloud density, (n w ), and are ordered by decreasing jet 
power from 10 46 erg s _1 down to 10 43 erg s _1 . All of these 
jets can disperse cloudy material for 77 > 77 cr it ~ 10 -3 5 
but for 77 < 10 -4 the jets would not strongly influence 
subsequent star formation. 

Next consider the BDG sequence. These simulations 
have a higher pisM/k and the same ratio of cloud to ISM 
densities and, thus, the same filling factor. The clouds 
are denser and harder to move, and in this sequence the 
critical value of the jet power to Eddington ratio, 77, is 



higher at around 10 -3 . 

We have also conducted simulations, in which the fill- 
ing factor is reduced by reducing the average density of 
the clouds for a fixed piswi/k. These are the sequences 
C f E f F f and B'D'G' plotted in Fig. [5] At high jet powers, 
these lower filling factor simulations all exhibit higher 
values of 77 cr i t , that is, it is more difficult for the jet to dis- 
perse the clouds when the clouds are more porous. This 
is the result of the non-thermal plasma being able to es- 
cape more readily through a more porous medium. On 
the other hand, the simulations F' and G' at jet powers 
of 10 44 ergs _1 occupy similar positions on the (v r ^ w )- 
Pjet diagram as their higher filling factor counterparts. 
In these cases the lower power jets do not break through 
as readily, and the confinement time of the non-thermal 
gas is sufficient to build up a large pressure, which ac- 
celerates the clouds. The results for even lower filling 
factors at all jet powers will be of interest but requires 
higher resolution simulations to prevent the clouds from 
becoming too pixelated. 

Another aspect of jet feedback, which is readily in- 
terpreted using this diagram is the relative importance 
for feedback of jets of different power. Compare the 
points H and C in Fig. [5] These points have approx- 
imately the same 77 cr it. However, the maximum value 
of the velocity dispersion, for which the jet in simula- 
tion H would satisfy the feedback criterion, (v ryW ) > a, 
is about 300kms _1 whereas the corresponding value for 
C is about 2, 300. Thus, we can make the qualitative 
remark that low-powered (Pj et < 10 43 ergs _1 ) jets are 
relevant to feedback in less massive hosts, that is, galax- 
ies in the early stages of hierarchical merging. 

3.3. Star formation in the over-pressured lobe 

Do stars form as a result of the excess pressure in the 
radio lobes? Let po = 10 _9 £>o,-9dyncm -2 be the ambient 
pressure surrounding a sphere with central number den- 
sity n c and temperature 10 4 T4 K. If magnetic fields are 
neglected, considerat ions of star formation focu s on the 
Bonnor-Ebert mass ( |Ebert 1955 Bonnor 1956) Mbe = 

M(T) or the 



2.0xl0 6 T 2 



6^2-1/2 ■ 



1.1S (kT/ flu) 2 G' 

50% larger critical mass M crit w 3.0 x 10 6 T 4 2 p~ ^ M de- 
rived from th e virial theore m for a uniform density gravi- 
tating mass (McCrea 1957). For a mass in excess of M cr it, 
say, we expect that pressure dr iven gravitational instabil- 
ity will lead to star formation (|Krumholz fc McKee) 2005 



Antonuccio-Delogu fc Silk||2008p . In the pre-jet phase o: 
our highest pressure simulations (p/k = 10 7 Kern" 3 ) the 
critical mass is approximately 2.5 x 10 6 M compared 
to the mass, 10 6 M , of the largest clouds with radius, 
R c ~ 25 pc. Hence, at this stage, the mass of the largest 
clouds is less than the critical mass, although not by a 
great amount. What is the effect on star formation when 
the ambient pressure increases as a result of the forma- 
tion of the high-pressured radio lobe? For example, in 
simulation F', the value of p/k, within the lobe rises from 
10 T K cm" ~ 3 to about 3 x 10 10 K cm , the critical mass 
decreases to about 5 x 10 4 M and jet induced starfor- 
mation may be feasible. 

However, the overpressured radio lobe is a highly tur- 
bulent environment and it is also helpful to compare 
the collapse timescale of the cloud, t c ~ (47rGp c ) _1//2 « 
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Figure 5. Maximum mean radial velocity of clouds against jet power for the simulations B-H of Table ^ The loci of constant 77, the ratio 
of jet power to Eddington luminosity are superposed. 



2 x 10 6 (n c /300 cm 3 ) x l 2 yr to the ablation time scale 
* ab i = 2^-^ * 8 x 10 4 (i? c /25pc)(^ abl /600kms- 1 )- 1 , 
where is the ablation velocity of cloud material. 
Clouds are destroyed on a time scales of a few xt ab i. If we 
allow the cloud density to increase by a maximal factor 
of 100 as a result of radiative shocks and adopt fiducial 
values for the other parameters we have t c ~ 2 x 10 5 yr 
and t a bi ~ 8 x 10 4 yr. For this extreme case one may ex- 
pect star formation to occur. On the other hand, a factor 
of 100 increase in density is high so that we expect that 
most clouds would be dispersed before they could grav- 
itationally coll apse to form stars. A sim i lar c onclusion 
was reached by |Antonuccio-Delogu fc Silk| ( |2008| ) in their 
study of the effect ol radio- lobes driven by jets with pow- 
ers ranging from 4 x 10 40 to 10 46 ergs -1 interacting with 
a single cloud with a radius of 10 pc. 

Given these estimates, under what circumstances 
would we expect jet-induced rather that jet-inhibited star 
formation? First, if the temperature of the clouds were 



lower than 10 4 K the critical mass for gravitational insta- 
bility (ex T 2 ) would also be lower. Second, we can also 
speculate that the clouds towards the edge of the radio 
lobe would not be so strongly affected by the lower level 
of turbulence in that region but they are still initially 
affected by the overpressure of the radio lobe bow-shock. 
Therefore it is possible that there may be some star for- 
mation in these regions. These issues are beyond the 
scope of this paper but they should be of interest in fu- 
ture work. 

3.4. Efficiency of kinetic energy transfer 

Panels a and b in Fig. [6] show the kinetic energy im- 
parted to the clouds in each simulation and the efficiency 
of the transfer of energy from the jet to the clouds, re- 
spectively. In each simulation the efficiency is quite high 
and, except for runs C and E', the efficiency is still in- 
creasing by the end of the simulation. All curves show 
a slight dip or inflection before steeply rising; this fea- 
ture occurs at the point where the jet plasma breaks 
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Figure 6. Evolution of various quantities in simulations B-H as functions of time. In all panels, the blue and green curves refer to 
simulations, for which = 1.0 and = 0.1, respectively, the solid curves refer to a volume filling factor fy = 0.42, and the dashed 
curves to fy = 0.13. (a) Kinetic energy of clouds. The dashed lines indicate the energy delivered by the jets with the jet power indicated 
on each curve, (b) Efficiency of energy transfer defined as the ratio of kinetic energy of clouds to the total energy delivered by the jet. (c) 
Density-averaged radial velocity of clouds, (d) Mechanical advantage, which is the ratio of the total momentum delivered by the jet to the 
total radial momentum in clouds. 



out of the confining clouds. (In the simulations of lower 
powered jets, the break out is smoothed out because the 
non-thermal plasma breaks out everywhere almost simul- 
taneously.) 

These curves show a decrease in efficiency with de- 
creasing value of the ratio of jet power to cloud density 
(Pjet/ (n w )). From the curves of (v r ^ w ) as a function of 
time shown in panel c of Fig. [6| we see that the maxima 
in (v r ^ w ) occur just before jet break-out, which releases 
some of the pressure in the quasi-spherical bubble. Com- 
pare the similarity of the evolution of E^[ njW / Pj et in B 
and E, D and F, G and H, B' and E' and D' and F'. 
These pairs of runs have the same values of Pj e t/ (n w ) 
and exhibit comparable maximum radial velocity disper- 
sion. 

3.5. Momentum transfer 

The rate of momentum transfer from the jet to the 
clouds can be an issue in this field, to which Krause fc] 
Gaibler (2010) drew attention. For exa mple, the range of 
cloud and jet parameters presented in Holt et al. (2006) 
indicates that the cloud momentum flux may exceed that 
of the jet. Hence, it is useful to examine the momentum 
budget to see what can be expected from the observa- 
tions. Momentum is not as straightforward as energy 
since the momentum in a volume surrounding a moving 
cloud is affected by the pressure (and magnetic stresses) 



integrated over the bounding surface of that volume. In 
this case the shocked jet gas provides a large pressure as 
a result of the energy flux in the jet. Therefore, in panel 
d of Fig. [6j we present the mechanical advantage (the 
ratio of the momentum in clouds to the jet momentum) 
as a function of time in each simulation. In each case 
the mechanical advantage exceeds unity. Note that the 
greatest mechanical advantage is achieved by the lowest 
power jets. Low powered jets "punch through" the warm 
clouds more slowly, and the shocked jet gas spreads more 
effectively through the volume of warm clouds and has a 
greater relative effect on the cloud momentum. 

3.6. Late-time evolution 

In each of the simulations, the jet plasma, after a suf- 
ficiently long time, escapes the cloud region. The more 
powerful the jet is, and the lighter the clouds are, the 
more centrally collimated the outbreak of the plasma is. 
We end the simulations when the jet plasma reaches the 
domain boundaries, and boundary effects become notice- 
able. We may, nonetheless, speculate on the subsequent 
evolution of the system. 

In the highest power jet simulations, we see the jet 
breaking through along its principal propagation axis 
(x = 0) and piercing the energy bubble that was inflated 
during the flood-and-channel phase of evolution. Except 
for runs G and H, the jets in all our simulations will 
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eventually emerge in a similar fashion, inflating a sec- 
ond, more classically shaped radio lobe beyond the first 
quasi-spherical bubble to a distance depending on the jet 
power. In runs G and H, the jet may be trapped indef- 
initely in the cloud region. While the common notion 
is that hierarchical radio lobe structures are prod uced 
by recurrent jet activity (Saikia & Jamrozy 2009), we 
note that they may also be produced by the random- 
ized emergence of radio plasma from an inhomogeneous, 
porous ISM. 

After the jet escapes along its principal propagation 
axis, the efficiencies in transfer of momentum and en- 
ergy to the clouds drop quickly. This stage is already 
seen in runs C and E', and will occur at later times in 
the other runs (except for G and H, if the jets in those 
runs remains trapped for the duration of their activity). 
While the clouds in our simulations are dispersed to high 
velocities, the host galaxy will retain most of the cloud 
material, except for a diffuse component comprising the 
high velocity tail in the distribution, which is concen- 
trated along the jet axis near the region where the cen- 
tral jet breaks out. By the end of the runs, the clouds in 
all simulations have reached and passed their maximum 
value of (v TjW ). 

4. DISCUSSION 

The suite of simulations presented in this paper 
strongly reinforce the importance of inhomogeneit y in the 
consideration of jet interactions with th e ISM (|Saxton| 
|et al.||2005l |Sutherland & Bicknelll 



1 



2007). Inhomogene- 



ity has several important consequences: (1) It affects the 
early morphology of the radio source as a result of the in- 
teraction of the jet and lobe with the obstructing clouds. 
(2) The radio source affects a much larger volume of the 
host galaxy because of the channeling of the jet flow in 
different directions. (3) The shocked jet clouds are left in 
th e wake of the non-thermal p lasma and, as first noted 
by Sutherland fc Bicknell| ( 2007 ) , would continue to emit 
shock excited line emission with the shocks driven by the 
high pressure gas in the non-thermal cocoon; this shock- 
excited emission is in addition to the emission that may 
be driven by photoionization by the nucleus. (4) The 
porosity of dense gas determines the ease with which a jet 
in a given host can disperse this gas, which is a potential 
source of new stars. In powerful sources, higher porosity 
gas is less easily dispersed but this trend is not evident for 
lower-powered (~ 10 44 ergs _1 ) jets. (5) Jets of all pow- 
ers can exert a considerable feedback effect on their host 
galaxies, although lower-powered jets only play a role in 
the lower velocity dispersion hosts. Brighter galaxies re- 
quire more powerful jets to disperse dense clouds. (6) 
The efficiency of transfer of kinetic energy from the jet 
to the dense gas is high. (7) The efficiency of transfer 
of momentum to the clouds is also high with mechanical 
advantages considerably exceeding unity. 

Inhomogeneity is therefore crucial when considering 
AGN feedback on the kiloparsec scale, both for the inter- 
pretation of radio and optical emission-line morphology 
in radio sources, which may be generating feedback, and 
for incorporating the effect of jet-mediated feedback on 
host galaxies of different size into large scale simulations, 
in which the resolution > 1 kpc. 

An important conclusion from these simulations is that 
jets with Eddington efficiency r] < 10 -4 are unlikely 



to have an effect on evolving galaxies when the pre- 
starforming gas exists in the form of clouds, which are 
relatively dense and cool compared to the hot ISM. This 
critical value of rj is relevant for clouds with a high fill- 
ing factor of 0.42 and a value of p/k — 10 6 for the ISM. 
We have shown that when the filling factor decreases or 
the pressure of the ISM increases the critical value of r] 
increases, and values of 7^ cr i t ~ 10 ~ 3 — 10 ~ 2 are not un- 
realistic. The precise values of ?7 cr it will have to await 
further higher resolution simulations with high porosity 
dense gas. 
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We are grateful to Prof. Mitchell Begelman for detailed 
comments on an initial version of the manuscript and 
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this work. We are grateful to our anonymous referee for 
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